STAR FORMATION HISTORY SINCE Z = 1.5 AS INFERRED 
FROM REST-FRAME ULTRAVIOLET LUMINOSITY 

DENSITY EVOLUTION 
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ABSTRACT 

We investigate the evolution of the universal rest-frame ultraviolet luminosity density from 
z = 1.5 to the present. We analyze an extensive sample of multicolor data {U'ab, Bab, Vab — 
24.5) plus spectroscopic redshifts from the Hawaii Survey Fields and the Hubble Deep Field. Our 
multicolor data allow us to select our sample in the rest-frame ultraviolet (2500 A) over the entire 
redshift range to z = 1.5. We conclude that the evolution in the luminosity density is a function 
of the form (1 + z)i-7±i o foj. ^ flat lambda {n^o = 0.3, Qxa = 0.7) cosmology and (1 + z)2 ''±i o 
for an Einstein-de Sitter cosmology. 



Subject headings: cosmology: observations — galaxies: distances and redshifts 
galaxies: formation — galaxies: luminosity function, mass function 



galaxies: evolution 



1. INTRODUCTION 

A major goal of observational cosmology is to 
understand the star formation history of the Uni- 
verse from the earliest epoch of structure forma- 
tion to the present. Much recent attention has 
focused on determining the contribution to the 
global history from the most distant sources; how- 
ever, the star formation history even at modest 
redshifts (z < 1) is not well determined and has 
recently undergone a revision. 

Early work by Madau et al. (1996) (later up- 
dated by Madau, Pozzetti, & Dickinson 1998) sug- 
gested that the global star formation as seen in the 
optical and ultraviolet (UV) had a strong peak 
around z = I and then fell very steeply at lower 
redshifts. The z < 1 data used in the analysis 
were taken from a paper by Lilly et al. (1996), 
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who used rest-frame near-UV luminosities derived 
from the /-selected Canada-France Redshift Sur- 
vey (CFRS, Lilly et al. 1995) to determine the 
comoving UV luminosity density from z = 1 to 
the present. These authors found the evolution to 
be a steep function of the form (1 + z)'*. How- 
ever, when Treyer et al. (1998) presented the 
first UV-selected constraints on the local inte- 
grated luminosity density, they found that their 
result was well above the optically-derived esti- 
mates. Sullivan et al. (2000) subsequently tripled 
the Treyer et al. UV-selected sample and con- 
firmed the higher local volume-averaged star for- 
mation rate. 

Cowie, Songaila, & Barger (1999, hereafter 
CSB) decided to reinvestigate the rest-frame UV 
luminosity density evolution to z = 1 using a large, 
extremely deep, and highly complete spectroscopic 
galaxy redshift survey. Their data enabled them 
to select objects based on the rest-frame UV mag- 
nitudes at all redshifts. The evolution found by 
these authors was a much shallower function of 
the form (1 -t- z)^-^. CSB suggested that the dif- 
ferences between their results and those of Lilly 
et al. could be accounted for by the 7-band se- 
lection of the Lilly et al. sample, which required 
large extrapolations to obtain UV colors, and by 
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the CFRS data not being deep enough to probe 
the flat segments of the luminosity function (LF), 
which meant that at redshifts near z = 1, reliable 
extrapolations to total luminosity density could 
not be made. 

In this paper, we expand on the work of CSB 
to more thoroughly investigate the rest-frame UV 
luminosity density evolution from 2: = 1.5 to the 
present. Our new galaxy sample is nearly twice as 
large as that used by CSB, and we explore various 
methods for constructing the UV LPs. The out- 
line of the paper is as follows. In § 2 we present 
our data sample and strategy and we investigate 
the U' number counts and redshift distribution. 
We also explore the redshift-magnitude relation- 
ship for the U', B, and V passbands. In § 3 we 
describe how we construct rest-frame UV LPs as 
a function of redshift from the U', B, and V data. 
In § 4 we utilize these LFs to infer the evolution 
of the global UV luminosity density with redshift. 
In § 5 we summarize our conclusions. Initially, we 
assume a flat lambda (f^mo = 0.3, Qxo = 0.7) cos- 
mology with Ho = 100 h km s~^ Mpc~^. Subse- 
quently, we investigate the effect of an Einstein-de 
Sitter (f2mO = 1-0, O^o = 0.0) cosmology on our 
results. 

2. OBSERVATIONS 

We analyzed a three passband subset {U' 
(3400 ± 150A), B, V) of an ongoing eight pass- 
band ([/', B, V, R, /, Z, J, HK') Hawaii imaging 
survey of four 6' x 2.5' areas crossing the Hawaii 
Survey Fields SSA13, SSA17, and SSA22 (Lilly, 
Cowie, & Gardner 1991) and the Hubble Deep 
Field (HDF) (WilUams et al. 1996). The B and 
V images were obtained using the Low-Resolution 
Imaging Spectrograph (LRIS; Oke et al. 1995) 
on the Keck 10 m telescopes and the UH8K CCD 
Mosaic Camera (Luppino 1997) on the Canada- 
France-Hawaii 3.6 m telescope. The U' data were 
taken with the ORBIT CCD on the University 
of Hawaii 2.2 m telescope. All magnitudes were 
measured in 3" diameter apertures and corrected 
to total magnitudes following the procedures de- 
scribed in Cowie et al. (1994). 

The sources analyzed in this paper contain and 
extend the CSB catalogs. The total number of 
galaxies with spectroscopic redshifts in the present 
U', B, and V samples are 403, 414, and 518, all 



to a survey limit magnitude (AB) of 24.5. These 
numbers can be compared to those in CSB (218, 
350, and 259 respectively) . Thus, our current 
datasct contains approximately double the num- 
ber of objects in the CSB U' and V samples. 

The great advantage of multiband data is that 
one can, for all redshifts, select galaxies based on 
their rest-frame UV magnitudes, thereby avoiding 
the uncertainties associated with selecting galaxies 
at longer (redder) wavelengths and then extrapo- 
lating to obtain their UV magnitudes. In addi- 
tion, the depth of the current dataset allows us to 
construct LFs to sufficiently faint absolute magni- 
tudes to constrain the faint-end slope. Moreover, 
by selecting galaxies based on their rest-frame UV 
magnitudes, we expect the relative shape of the in- 
ferred UV LFs to be minimally sensitive to the ef- 
fect of interstellar dust; i.e., we expect the relative 
forms of our derived LFs (and hence the relative 
values of our luminosity densities) to be subject 
only to evolution in the amount or properties of 
galactic dust with redshift and not to its effect in 
absolute terms. (Note that any extinction due to 
dust in the source galaxy should not be confused 
with the tiny extinction corrections due to inter- 
stellar dust) 

As discussed in CSB, one is free to choose the 
rest-frame wavelength at which to compute LFs 
and the UV light density. For our dataset, 2500 A 
provides a sensible compromise between our wide 
range of redshifts and large number of galaxies in 
each passband. 

Because there are missing or unidentified galax- 
ies in each sample, we construct LFs by either 
omitting them (i.e., we assume they are imidcn- 
tified because they are unusual in some way, e.g., 
at very high redshift), or by assuming they are 
distributed in redshift in exactly the same manner 
as the identified galaxies (i.e., we assume they are 
identical to the measured galaxies but have been 
missed for some trivial reason). It is likely that 
the true LF lies between these two possibilities. 
In subsequent sections we refer to the former case 
as minimal and the latter case as incompleteness- 
corrected. The [/', B, and V samples are 88%, 
90%, and 83% complete, respectively, so the cor- 
rection is not a huge factor in any case. 
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2.1. The U' Sample 

Figure 1 shows number counts versus appar- 
ent magnitude for our U' sample. The symbols 
indicate the number counts in each of the four 
fields. We assume a 1 cr Poisson uncertainty for 
each field. The solid line shows the best-fit to the 
counts for galaxies with magnitudes between U' of 
22.0 and 24.0 (slope of 0.61 ±0.06 and intercept at 
— 10.34 ± 1.49). The uncertainties were calculated 
from the field- to-field variations. 

A number of groups (e.g., Williams et al. 1996; 
Pozzetti et al. 1998; Gardner, Brown, & Ferguson 
2000) have measured deep galaxy counts at 3000 A 
from the HST Hubble Deep Field imaging survey, 
but their counts are somewhat deeper than those 
presented here, so it is difficult to make a direct 
comparison. The sample most similar to ours was 
that obtained at the Palomar 5 m by Hogg et al. 
(1997). We overplot their data (ffiled diamonds) 
on Fig. 1 for comparison, after converting their 
magnitudes to U'ab by adding 0.79 mag. The 
agreement is generally good. 

In Fig. 2 we show the total number of objects 

in our U' sample versus redshift for two apparent 
magnitude bins. The upper panel shows the num- 
ber of galaxies versus redshift for all galaxies with 
U' between 22.5 and 23.5, and the lower panel is 
for galaxies with U' between 23.5 and 24.5. 

Wilson et al. (2001) found that a good model 
for the redshift distribution (at least at I and V) 
is provided by 



p{z) = 0.5z'^ exp{—z/ zq)/ zq'^ 



(1) 



where p(z) x dz is the probability of finding a 
galaxy in the redshift interval z + dz (the mean 
redshift is z = Szq and the median redshift is 
-^median = 2.67zo). A nicc property of (1) is that 
there is only one free parameter, the redshift scale 
parameter, zq. The solid lines overlaid on Fig. 2 
show the best-fits to the model {zq = 0.19 and 
zo = 0.28), normalized to the total number of 
galaxies in each sample. 

2.2. Magnitude-Redshifb Dependence 

In this section we investigate the dependence of 
median redshift on apparent magnitude and wave- 
length. In Fig. 3 we show redshift versus magni- 
tude for the U' sample. The symbols denote the 
field in which the galaxy was observed. Note that 




Fig. 1. — Number counts versus apparent mag- 
nitude for the U' sample. The four fields are as 
indicated by the key. The solid line shows the 
best-fit to the counts between U' of 22.0 and 24.0 
(slope of 0.61 and intercept at —10.34). The filled 
diamonds are the counts from Hogg et al. (1997). 




Redshift 

Fig. 2. — Redshift distributions for the U' sam- 
ple. The upper panel shows N{z) for a one mag- 
nitude wide band (22.5 < U' < 23.5) in appar- 
ent magnitude. At this depth the data are 97% 
complete. The solid line shows the best-fit to the 
model (Eq. 1) with redshift scale parameter 0.19. 
The lower panel shows N{z) for 23.5 < U' < 24.5. 
At this depth the data are 87% complete and best- 
fit with redshift scale parameter 0.28. 
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Table 1 

Median redshifts and galaxy numbers as a function of apparent magnitude and 
wavelength (the uncertainties are 1 fj poissonian limits). 
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Fig. 3. Galaxy rcdshift versus apparent magni- 
tude for the U' sample. The key is as in Fig. 1. 
The solid line shows the median redshift and the 
1 (7 Poisson uncertainties (Table 1). 



the SSA13, SSA22, and HDF fields are - 90% 
complete to an AB limiting magnitude of 24.5, and 
the SSA17 field is similarly complete to 23.5. 

Table 1 quantifies Fig. 3, giving the median rcd- 
shift with ± 1 a Poisson uncertainties (Gehrels 
1986) and the number of objects in each half- 
magnitude interval as a function of apparent mag- 
nitude. We note that the median redshift obtained 
from the parameterized fit to the data (Eq. 1) for 
the 22.5 <U' < 23.5 interval is 0.51, and for the 
23.5 <U' < 24.5 interval, 0.75. The values are in 
good agreement with those calculated directly in 
Table 1. 

In Table 1 we also calculate the median red- 
shift for the B and V samples. Median redshift as 
a function of apparent magnitude can potentially 
be used to constrain galaxy evolution models, e.g. 
if we compare with the B-band predictions from 
the merger model proposed by Carlberg (1992, his 
Table 2), we find that his median redshift values of 
0.44 for B = 23 and 0.55 for B = 24 are lower than 
our values, suggesting that more pure luminosity 
evolution might have occurred than was proposed 
in his model. 

3. THE REST-FRAME ULTRAVIOLET 
LUMINOSITY FUNCTION 

For any given galaxy, i, at redshift, z, the equa- 
tion relating the absolute and apparent magni- 
tudes is given by 



Mf 00 = mi-5logdL{zi)-25+2.51og{l+z)+dK{z^) 

(2) 
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where m is the observed magnitude at the red- 
shifted wavelength and dL{z) is the luminosity dis- 
tance in Mpc. dK{z) is given by 



dK{z) = 2.5 log 



10 



./;,2o()(.)A X (1 + ^)] 
/.[2500Ax (l + 0c)] 



(3) 



where fv is the spectral energy distribution of the 
galaxy and Zc is the redshift corresponding to the 
center of the band. dK{z) is a differential K- 
correction to account for each sample containing a 
range of redshifts and hence a range of rest-frame 
wavelengths around 2500 A. It is generally small 
and is obtained by interpolation from the neigh- 
boring passbands. 

3.1. Luminosity Functions from the V^ax 
Method 

We used two methods to construct LFs: the 
traditional T^ax method described by Schmidt 
(1968), Felten (1976), and Ellis et al. (1996) and a 
new method recently suggested by Page & Carrera 
(2000). In the Fmax method, the number density 
of galaxies in the redshift range [zi , Z2] with abso- 
lute magnitude M is given by 



0(M) dM = J2 



1 



Knax(M) 



(4) 



where the sum is over the galaxies in the magni- 
tude interval M ± dM/2. Vmax(-M) is the max- 
imum total comoving volume within which each 
galaxy (as defined by its apparent magnitude and 
redshift) would remain detectable within survey 
limits. The uncertainty for each magnitude inter- 
val is conventionally calculated from 



E 



1 



(Kna.(M))2 



(5) 



(Marshall 1985; Boyle, Shanks, & Peterson 1988). 
This expression weights each observation by its 
contribution to the sum. However, it assumes 
Gaussian statistics, which is not ideal for bins at 
the bright or faint end of the LF where only a 
small number of objects contribute to the sum. 

Figure 4 shows the 2500 A rest-frame LF for 
our three redshift bins: z = 0.2 — 0.5, con- 
structed from the U' sample (upper panel), z = 



0.6 — 1.0, constructed from the B sample (cen- 
ter panel), and z = 1.0—1.5, constructed from 
the V sample (lower panel). At redshifts z = 
0.35 ± 0.15, 0.80 ± 0.20 and 1.25 ± 0.25, the U' 
(3400 A), B (4500 A) and V (5500 A) samples cor- 
respond to rest- wavelengths of 2519^253' 25OOI250 
and 2444^244 A. The open circles denote the min- 
imal function, and the solid circles denote the 
incompleteness-corrected function. The number of 
galaxies used in the construction of each LF was 
121 ([/'), 119 (B), and 59 (V). 

3.2. Luminosity Functions from the PC 

Method 

The Page & Carrera (2000) method is very sim- 
ilar to the V^ax(-^) method in that it results in a 
binned differential LF, but the advantage is that 
it more accurately determines the LF at the faint 
end. The maximum redshift at which any galaxy 
may be found is a constantly varying function de- 
termined by the flux limit of the survey. The 
14nax(Af) method assumes that the redshift is a 
constant for any given absolute magnitude bin. 
However, by dividing each magnitude bin into a 
series of steps, calculating Knax(Af) for each in- 
terval, and then integrating over the magnitude 
bin, one can obtain a more precise estimate of the 
maximum volume to which each object could be 
observed in the survey. Estimates of (f){M) ob- 
tained by either of the two methods should agree 
at the bright end of the LF where objects are much 
brighter than the flux limit. The two estimates 
should also agree in the case of very fine magnitude 
intervals, such that the widths of the magnitude 
bins tend to zero. 

We denote the LFs obtained using the Vpc 
method by diamonds in Fig. 4. The open triangles 
show the minimal function and the solid triangles 
show the incompleteness-corrected function. The 
diamonds have been offset from the circles of the 
14iax method by 0.15 mag for clarity. The ± 1 
a Poisson uncertainties (Gehrels 1986) in the LFs 
are appropriate for small numbers of objects per 
magnitude interval. The LFs obtained from the 
two methods are very similar. 

Both Vmax(-M) and Vpc{M) give unbiased es- 
timates of 4>{M) only if galaxy clustering can be 
neglected. It is easy to imagine how a nearby ex- 
cess of clustering could bias the estimator: such 
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-22 -20 -18 -16 -U 



M(2500) + 5 log(h) 

Fig. 4. — The 2500 A rest-frame luminosity fimc- 
tion for three redshift bins: z = 0.2 — 0.5, from the 
V sample (upper panel), z = 0.6 — 1.0, from the 
B sample (center panel), and z = 1.0 — 1.5, from 
the V sample (lower panel). An fimo = 0.3, r^Ao = 
0.7 cosmology is assumed. The open circles de- 
note the minimal function, and the solid circles 
denote the incompleteness-corrected function ob- 
tained using the Knax method. The diamonds 
(offset by 0.15 mag for clarity) show the minimal 
(open) and incompleteness-corrected (solid) func- 
tions obtained using the Vpc method. See text 
for details and an explanation of the uncertainties. 
For this dataset the luminosity functions obtained 
from the two methods are very similar. 



an excess of intrinsically faint galaxies would cause 
the LF to be too steep at the faint end. Although 
the effects of clustering are mostly of concern in 
pencil-beam surveys, wc tested two maximum like- 
lihood alternatives, which should be less sensi- 
tive to clustering, on our U' sample. These were 
the Schechter fit estimator suggested by Sandage, 
Tammann, & Yahil (1979) and the stepwise es- 
timator of Efstathiou, Ellis, & Peterson (1988). 
In all cases, the LFs were very similar to those 
obtained using Vm^x, giving us confidence in the 
robustness of our results. We also note that the 
LFs in this paper were obtained using code written 
independently from that used in CSB. 

3.3. Schechter Parameterization 

Wc then assumed that each LF could be pa- 
rameterized by a Schechter function 

^{M) dM = fc0*e'=(«+i)(^*-^)e-«''"*""* dM 

(6) 

where A; = fin 10 (Schechter 1976). We solved 
for the best-fit Schechter parameters assuming two 
fixed faint-end slopes that likely bound the range 
of faint-end slopes : a = —1.0 and a = —1.5. Ta- 
ble 2 shows the best-fit absolute magnitude at the 
knee, M*, and normalization, (j)* for these values 
of fixed faint-end slope. In Fig. 5 the open circles 
again show the minimal function; the solid circles 
the incompleteness-corrected function assuming a 
flat lambda cosmology as in Figure 4. Overlaid on 
Fig. 5 arc the best-fitting Schechter functions as- 
suming a = —1.0 or a = —1.5. The solid portion 
of the line shows the magnitude range utilized in 
the fit. 

3.4. Effect of Cosmology 

We then investigated how one's choice of cos- 
mology affects the LFs. We reconstructed LFs 
again using the V^ax(-M) method but this time 
assuming an Einstein-de Sitter cosmology. For 
comparison, in Fig. 5, the open triangles show 
the minimal function; the solid triangles the 
incompleteness-corrected function assuming this 
cosmology. As in Figure 4, the triangles have 
been offset by 0.15 mag for clarity. 

Table 3 shows the best-fit Schechter values of 
absolute magnitude at the knee (M*) and normal- 
ization {(/)*) for the Einstein-de Sitter cosmology 
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Table 2 

schechter function parameter fits for minimal and incompleteness-corrected points for 
flat lambda cosmology. the columns show the best-fit values of absolute magnitude at 
the knee (m*), normalization {^*), and reduced for fixed faint- end slope (a) of -1.0 

AND -1.5. 



a = -1.0 a = -1.5 



Sample 


Corrected? 


M* 






M* 






Range 


U' 


No 


-18.17 


0.0142 


1.90 


-19.86 


0.0024 


2.51 


[-19.25 


-14.75] 


U' 


Yes 


-18.08 


0.0161 


1.25 


-19.64 


0.0031 


1.58 


[-19.25 


-14.75] 


B 


No 


-18.46 


0.0178 


0.55 


-18.96 


0.0089 


0.86 


[-21.25 


-17.25] 


B 


Yes 


-18.34 


0.0232 


0.20 


-18.81 


0.0123 


0.36 


[-21.25 


-17.25] 


V 


No 


-18.30 


0.0212 


0.87 


-18.53 


0.0183 


0.75 


[-20.75 


-18.25] 


V 


Yes 


-18.12 


0.0348 


0.69 


-18.35 


0.0306 


0.58 


[-20.75 


-18.25] 



Table 3 

schechter function parameter fits for minimal and incompleteness-corrected points for 
EiNSTEiN-DE Sitter cosmology. The columns show the best-fit values of absolute 

MAGNITUDE AT THE KNEE (M*), NORMALIZATION ((/)*) AND REDUCED FOR FIXED FAINT- END SLOPE 

(a) OF -1.0 AND -1.5. 



a = -1.0 a = -1.5 



Sample 


Corrected? 












X^l^ 


Range 


U' 


No 


-17.64 


0.0322 


0.41 


-19.02 


0.0070 


0.88 


[-18.75 


-14.75] 


U' 


Yes 


-17.52 


0.0378 


0.22 


-18.75 


0.0094 


0.41 


[-18.75 


-14.75] 


B 


No 


-17.86 


0.0512 


0.26 


-18.23 


0.0318 


0.25 


[-20.75 


-17.25] 


B 


Yes 


-17.75 


0.0643 


0.09 


-18.13 


0.0405 


0.05 


[-20.75 


-17.25] 


V 


No 


-17.63 


0.0812 


0.19 


-17.85 


0.0699 


0.21 


[-20.25 


-17.75] 


V 


Yes 


-17.54 


0.1100 


0.04 


-17.78 


0.0932 


0.04 


[-20.25 


-17.75] 
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(again assuming fixed faint-end slopes of —1.0 and 
— 1.5). From Figure 5 and from Tables 2 and 3, 
it is clear that the larger distance/volumes associ- 
ated with a cosmological constant cause M* and 
(j)* to decrease compared to the best-fit parameters 
in the Einstein-de Sitter case. 



4. REST-FRAME UV LUMINOSITY 
DENSITY EVOLUTION 

The LFs in § 3 can now be used to calculate the 
rest-frame ultraviolet luminosity densities, C, with 
redshift. One approach is to choose a magnitude 
limit and to sum the LF over the magnitude bins 
directly using 

-|^Q-0.4M 

>CDirect = 4.4x1020^ —h erg s-1 Mpc 

(7) 

An alternative approach, which we adopt, is to 
choose a faint-end slope (we use either a = —1.0 
or a = —1.5) and to integrate the LF analytically 
using the best-fit Schechter parameters from Ta- 
bles 2 and 3 

poo 

-Cschechter= / L(j,{L)dL = L^(l)^T{a + 2) (8) 

Jo 

giving 




-22 -20 -18 -16 -14 



M(2500) + 5 log(h) 



/:Schechter = 4.4 X lO^" X lO'^-^ ^*(/.^r(a + 2) h 

(9) 

in units of ergs s~^ }iz~^ Mpc"'' . Although this 
method involves integrating over all luminosities, 
fainter galaxies have a rapidly decreasing contri- 
bution to the total luminosity density, and thus 
the two methods give similar results. We cal- 
culate the luminosity density for both the mini- 
mal and incompleteness-corrected cases assuming 
firstly a faint-end slope of a = —1.0 and then 
a = —1.5. The resulting luminosity densities for 
the flat lambda and Einstein de-Sitter cosmologies 
are shown in Table 4. 

Luminosity density evolution with redshift is of- 
ten parameterized as a power-law, £ oc (1 -|- z)^. 
In Fig. 6, we show log^Q (luminosity density) versus 
logiQ(l -|- z) using the values from Table 4. As in 
Fig. 5, we use circles to denote the the flat lambda 



Fig. 5. — Circles as in Fig. 4 for a Q,mo = 
0.3, rixo = 0.7 cosmology. Overlaid are the best- 
fitting Schechter functions assuming either a = 
— 1.0 or a = —1.5. The solid portion of the 
line shows the magnitude range utilized in the 
fit. Also shown and denoted by the triangles 
(offset by 0.15 mag for clarity) are the minimal 
(open) and incompleteness-corrected (solid) func- 
tions obtained for an Einstein-de Sitter cosmology. 
Sec text for details. See Tables 2 and 3 for best 
Schechter function fit parameters. 
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Table 4 

CoMoviNG 2500A UV Luminosity Density.'' 



logio (Luminosity Density) 



Sample 


Redshift 


Corrected? 


Flat Lambda 
a = -1.0 a = -1.5 


Einstein-de Sitter 
a = -1.0 a = -1.5 


U' 


0.35 ±0.15 


No 


26.058 


26.211 


26.202 


26.340 


U' 


0.35 ±0.15 


Yes 


26.077 


26.234 


26.223 


26.360 


B 


0.80 ±0.20 


No 


26.272 


26.420 


26.491 


26.681 


B 


0.80 ±0.20 


Yes 


26.339 


26.500 


26.546 


26.746 


V 


1.35 ±0.25 


No 


26.284 


26.561 


26.600 


26.871 


V 


1.35 ±0.25 


No 


26.428 


26.712 


26.695 


26.968 



''Units are h ergs s ^ Hz Mpc 



cosmology and triangles to denote the Einstein-de 

Sitter cosmology. The open symbols denote the 
minimal case , and the solid symbols denote the 
incompleteness-corrected case. 

We solved for the best-fit power-law exponent, 
(3, in each case. We used the mean of the luminos- 
ity densities obtained assuming faint end slopes of 
a = — 1.0 and a = —1.5 as our best estimate, with 
the extreme values as estimates of the uncertainty. 
Table 5 gives the best-fit exponent and uncertainty 
function of completeness-correction and cos- 
mology. For the flat lambda cosmology we found 
a best-fit exponent of 1.44 ± 0.63 in the minimal 
case and a best-fit exponent of 1.95 ± 0.65 in the 
incompleteness-corrected case. For the Einstein 
de-Sitter cosmology we found a best-fit exponent 
of 2.22 ± 0.62 in the minimal case and a best- 
fit exponent of 2.54 ± 0.62 in the incompleteness- 
corrected case. Thus, depending on the choice of 
completeness correction, we conclude that lumi- 
nosity density evolves as (l+z)^'^^^ ° in the fimo = 
0.3, ^AO = 0.7 cosmology and as (1 + 2;)2-4±i o 
the Einstein de-Sitter cosmology. The two solid 
lines overlaid on Fig. 6 show the best-fit solu- 
tions for each cosmology in the incompleteness- 
corrected case. The Einstein de-Sitter value is 
slightly steeper than that obtained by CSB (1.3 
for a = —1.0 and 1.7 for a = —1.5) but consistent 
within the uncertainties. 

Finally, in Fig. 7 we compare our values of lu- 




o This Work (Q^ = 0.3, = 0.7) 
A This Work (0^ = 1.0. = 0.0) 



0.2 0.3 
log (l+z) 



Fig. 6. — log (Luminosity density) versus 
logiQ(l + z). The open circles denote the 
minimal case and the solid circles denote the 
incompleteness-corrected case, both for the flat 
lambda cosmology. The vertical bars show the un- 
certainties in luminosity density caused by assum- 
ing a = —1.0 (lower) or a = —1.5 (upper). The 
triangles show same and are for an Einstein de- 
Sitter cosmology. They have been offset slightly 
for clarity. Also shown arc the best-flt power-law 
to the incompleteness-corrected values for each 
cosmology. See § 4 for a discussion of the best 
power-law fits and Table 5 for the values. 
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Table 5 

Best-fit Power Law Exponent /? 



Corrected? 


Flat Lambda 


Einstein-de Sitter 


No 


1.44 ±0.63 


2.22 ±0.62 


Yes 


1.95 ±0.65 


2.54 ±0.62 



m (N 

m 



E in 

g CM 

% 



A This Work 

» Sullivan et al. 2000 

* Pascarelle et al. 1998 

* Connolly et al. 1997 
n Lilly et al. 1996 

I I I I I I I I L 



0.2 0.3 
log (l+z) 



0.4 



Fig. 7. — log (Luminosity density) versus 
logio(l + z) for an Einstein de-Sitter cosmology. 
The triangles are as in Fig. 6 and denote the min- 
imal (open) and incompleteness-corrected (solid) 
case. The vertical bars show the uncertainties in 
luminosity density caused by assuming a = —1.0 
or a = —1.5. For comparison we have included 
data points from Lilly et al. (1996, open squares), 
Connolly ct al. (1997, open star), Pascarelle ct al. 
(1998, open cross) and Sullivan et al. (2000, open 
diamond). 



minosity density to the values obtained by other 
surveys. Other groups have previously assumed an 
Einstein-dc Sitter cosmology and therefore should 
be compared to the triangles from Fig. 6. To 
convert the low-rcdshift value obtained by Sulli- 
van et al. (2000) (from the FOCA2000 balloon- 
born survey) we used their best-fit Schechter val- 
ues (their Table 3) and converted their magnitudes 
to AB magnitudes using a 2.29 mag offset; we also 
converted from a rest-frame of 2000 A to 2500 
A using a A^ "'^ power-law, as suggested by Fig. 4 
of CSB. As mentioned in § 1, the Sullivan et al. 
UV-selected sample results in a higher value of in- 
tegrated luminosity density for the local universe 
than previous optically-derived estimates. 

In comparing with the surveys of Pascarelle, 
Lanzetta, & Ferndndez-Soto (1998), Connolly 
et al. (1997), and Lilly et al. (1996), we again 
converted to a rest-frame of 2500 A using a 
A"'^ "'^ power-law. We also converted to Hq = 
100 h km s~^ Mpc~^, where necessary. Con- 
nolly et al. and Pascarelle et al. both calculate 
their luminosity density from the small Hubble 
Deep Field North-proper (HDF-N) using photo- 
metric redshift estimates. The results of Connolly 
et al. are more directly comparable to ours since 
they measure at a rest-frame of 2800 A whereas 
Pascarelle et al. measure at a rest-frame of 1500 
A. Connolly et al. assume a faint end slope of 
a = —1.3. Both obtain somewhat higher lumi- 
nosity densities than do we, although the values 
are consistent within the imcertainties. As dis- 
cussed at length in the Appendix to CSB, much of 
the discrepancy between our results and Connolly 
et al. and Pascarelle et al. may be attributable to 
the slightly higher number counts in the HDF-N 
versus other fields. 

Lilly et al. calculated their luminosity den- 
sity from the Canada- France Redshift Survey, and 
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found a best-fit exponent of /3 = 3.9 ± 0.75. As 
discussed in § 1, this is somewhat steeper then 
the value of /? = 2.5 ±1.0 that we obtained. Prom 
Fig. 7, wc conchidc that the steeper vahie obtained 
by LiUy et al. is most hkely due to a combination 
of their 2; ~ 1 luminsity density estimate being 
rather higher than ours, and their use of a low 
(optically-derived) estimate of the local luminsity 
density. 

In closing, it is important to add one caveat con- 
cerning the effect of interstellar dust on our conclu- 
sions. In this paper we assumed that any extinc- 
tion would suppress UV emission uniformly. This 
corresponds to applying a constant correction fac- 
tor to the LFs, and does not affect the luminosity 
density slope inferred from Fig. 7. Some authors 
have suggested that extinction may be luminosity 
dependent (Adelberger & Steidel 2000; Sullivan 
et al. 2000; Hopkins et al. 2001). If this is the 
case, the higher redshift luminosity functions con- 
taining greater contributions from brighter galax- 
ies would see larger extinction corrections, possi- 
bly flattening the slope observed in Fig. 7 from 
a steeper value. The satisfactory resolution of 
the complex role of dust and the validity of these 
claims will require further investigation with larger 
samples. 

5. SUMMARY 

We investigated the evolution of the universal 
rest-frame luminosity density from z = 1.5 to the 
present. The availability of both multicolor data 
and highly complete spectroscopy enabled us to 
select galaxies based on their rest-frame ultravi- 
olet color, minimizing potential sources of error 
such as large if -corrections and interstellar dust. 
Our large, deep sample allowed us to constrain 
the faint-end of the luminosity function with con- 
fidence, even at the highest redshift interval of 
z = 1.25 ± 0.25. Assuming analytic Schechter 
forms for our luminosity functions and using likely 
extremal faint-end slope values of a — —1.0 and 
a = —1.5, we constrained the relative luminosity 
density as a function of redshift. We concluded 
that, in an fimo = 0.3, 51 ao = 0.7 [Einstein-de Sit- 
ter] Universe, the evolution in the luminosity den- 
sity follows ^)i-7±i o [(1 + 2)2.4±1.0j yiopg 

from z = 1.5 to the present, implying that rather 
more star formation has occurred in recent times 



than was previously suggested. 
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